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Abstract 



The study for lepton flavor violation combined with the neutrino oscillation 
may provide more information about the lepton flavor structure of the grand 
unified theory. In this paper, we study two lepton flavor violation processes , 
t — ► ^7 and Z — > r/i, in the context of supersymmetric SO (10) grand unified 
models. We find the two processes are both of phenomenological interest. 
In particular the latter may be important in some supersymmetric parameter 
space where the former is suppressed. Thus, Z-dacay may offer another chance 
for looking for lepton flavor violation. 



I. INTRODUCTION 



Super-Kamiokande data on the atmospheric neutrino anomaly presents a strong evidence 
for the existence of neutrino oscillations. The anomaly can be explained by — v T oscillation 



with 5m| 3 = (1 — 8) x 10~ 3 eV 2 and a large mixing angle sin 2 (20 MT ) =0.8 — 1 |],|2|. In 
addition, the long standing solar neutrino deficit can also be interpreted as another type 
of neutrino oscillations ||[|]. Assuming that the LSND anomaly || will finally disappear 
all these observations about neutrinos can be accommodated in a model with three very 
light active left-handed neutrinos. Many such models have been proposed to explain the 
measured neutrino parameters since Super-Kamiokande first published the data ||. Among 
them a natural explanation of the neutrino masses is provided by grand unified models which 
interpret the very light neutrino masses compared with the quarks and charged leptons by 
the see-saw mechanism 0. Especially in SO(10) grand unified models, in which the right- 
handed (RH) neutrinos have masses of the order about unification scale and the lepton masses 
and the quark masses are related, correct neutrino masses can be obtained j5|. The measured 
neutrino masses are even regarded as a new evidence supporting the grand unification idea 
|§. As been argued in Ref. || a value for m{v T ) ~ falls into a natural range predicted 
by a grand unified model based on either a string-unified G(224) model or a SO(10) grand 
unified model. Furthermore, grand unified theory (GUT) has the advantage that it relates 
the neutrino problem, its masses and mixings, with charged lepton and quark masses and 
mixings into a large fermion flavor problem and thus gives more definite predictions. 

Conversely, the parameters measured in the neutrino sector can provide a window on 
the grand unified theory study and they also give important feedback on the problems of 
quarks and charged lepton masses as they are all related in the GUT. A direct inference of 
neutrinos being massive is the existence of a Kabayashi-Maskawa like matrix for the lepton 
sector. However, the processes violating lepton flavors due to this matrix is too small to 
be observed because of the tiny neutrino masses. In the supersymmetric GUTs the high 
energy lepton flavor violation (LFV) interactions may leave trace in the mass matrices of 
scalar partners of leptons by renormalization effects, which drives low energy LFV processes 
since the mass matrix of the charged leptons and that of sleptons can not be diagonalized 



2 



simultaneously. Therefore, the study of the LFV process may provide important information 
on the flavor structure of the supersymmetric GUTs. Various LFV processes have been 



studied in different theoretical frameworks by a number of authors [10-14 



There are two purposes of this work. One is to study the LFV processes in the context of 
supersymmetric SO(10) grand unified models with a "lopsided" texture for the mass matrices 
of the down quark and charged lepton which has been advocated by a number of authors to 
accommodate the large — u T mixing and small jn|[n].[| . If the experimental sensitivity 



on t — > /ry can reach down to 1 x 10~ 9 [|1(J, not only this process can be observed but it can 
even be used to discriminate different fermion textures in the grand unified models. Another 
purpose of this work is to present the analytic formulas for the branching ratio of the process 
Z — > Tfi in the minimal supersymmetric standard model (MSSM) and the numerical study 
in the class of models considered by us. This calculation seems absent in the literature. This 
calculation is triggered by the Giga-Z option of the Tesla project which may expect the 10 9 
Z bosons at resonance jD]]. The upper limit for the branching ratios could be improved 
down to BR(Z — > ^t^) < f x 2.2 x 1CT 8 with / = 0.2 ~ 1. We find this process is quite 
interesting phenomenologically since in some supersymmetric parameter space r — ► /ry being 
suppressed, this process can be important however. 

The paper is arranged as follows. In sec II we discuss the origin of the LFV interactions 
in a supersymmetric grand unified model and the Renormalization Group Equations(RGE). 
In sec III we present the formulas for the branching ratios of the two processes. In Sec IV we 
briefly introduce a grand unified model proposed in Ref [15] that we used in our calculations. 



The numerical results are presented in Sec V and we summarize and give conclusions in Sec 
VI. 
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II. RENORMALIZATION GROUP EQUATIONS AND LOW ENERGY 
SUPERSYMMETRIC SPECTRUM 



A. Origin of Lepton Flavor Violation 



If the SM is extended with massive and non-degenarate neutrinos, LFV processes may be 
induced. However, such processes are highly suppressed due to the smallness of the neutrino 
masses. The branching ratio is proportional to 5ml/ which is hopeless to be observed 
HI. When supersymmetry enters the theory the scene changes completely. The LFV may 



also be induced through the generation mixing of the supersymmetric soft breaking terms 
of the lepton sector. However, arbitrary mixing of these soft terms in the MSSM are not 
predictive. In our calculations we adopt the supergravity mediated supersymmetry breaking 
mechanism to produce universal soft terms at the GUT scale Mqut, because non-universal 
soft terms at M GUT may produce too large low energy LFV observable effects |T0| , pT| . The 
tree level universal soft terms may induce non-diagonal terms at low energy by radiative 
corrections including LFV interactions at high scale. Our procedure includes calculating the 
low energy supersymmetric soft terms which is not generation universal now by integrating 
the RGEs and then calculate the LFV branching ratios induced by these non-universal soft 
terms. 

For a supersymmetric SO(10) grand unified model, the structure below Mqut where the 
grand unification has been spontaneously broken is the same as the MSSM supplemented 
with MSSM singlet RH neutrino superfields. The superpotential of the lepton sector is now 

W = PjHiliNj + ffHiLiEj + -M^NiNj + ^HxH 2 (2.1) 

where f v and fi are the Yukawa coupling matrices, M is the RH neutrino mass matrix, i, 
j are the generation indices. Anti-symmetric tensor e ab is implict to contract the SU(2) 



doublets with e 12 = —1. In general, f u and fi can not be diagonalized simultaneously, which 
is the origin of LFV interactions. Diagonalize f u and // by bi-unitary rotations, 

ff = uIWr 

ft = Vtf v V R (2.2) 

where Ul,r, Vl,r are all unitary matrices. Then define 

V D = U{V L (2.3) 

which is analog to the KM matrix Vkm in the quark sector. Vd is crucial for LFV pro- 
cesses. The RH neutrino masses are not much lower in order compared to Mqut in 
SO (10) grand unified models. After the RH neutrinos are decoupled and H 2 gets VEV 
V2 of weak scale we get three light left-handed (LH) Majorana neutrinos with mass matrix 
m u = —(fvV2)M~ 1 (f u V2) T by see-saw mechanism. Suppose m v m\, is diagonalized by V£, 
then the matrix 

yMNS = jjfyu (2A) 

determines the neutrino oscillation parameters. In the published grand unified models which 
emphasize the neutrino oscillations the large — v T mixing angle in V MNS is mainly coming 
from W L 0. Thus, on the one hand, we may expect large rates for LFV processes due to 
large fi — t mixing in U^, which at the same time causes large — v T mixing in V MNS 
for neutrino oscillations, and, on the other hand, difference between Vd and V MN may be 
found by neutrino oscillations and LFV processes. If so, important information about GUT 
structure may be derived then. 



The soft breaking terms for the lepton sector is [19|, 



C soft = -m\H\H x - m 2 H2 H\H 2 - {m\) ir L\Lj - {m\fk*R 3 - {mlfv*^ 



+(A l/ f u )^H 2 L i D j + h.c.) 



(2.5) 



where i, j are generations indices. At M GUT we assume the universal conditions, 




m 2 Hl = ml 



(2.6) 



2 2 2 

m 6 = m r , = m n 



(2.7) 



A v — A\ 



o 



(2.8) 



Fig. 1 gives the explanation of the occurrence of low energy LFV processes. The non- 
diagonal Yakawa couplings induce z> M — v T mixing through loop effects. This high energy 
process can be running down by integrating the RGEs to low energy. Thus the study of LFV 
processes induced by non-universal soft terms can actually reveal high energy fermion flavor 
struture. In the basis where the /; is diagonal we can get the non-diaganol scalar mass in 
the first order of approximation as 



where in the diagonalized Yukawa matrix fl only the (3,3) element f U3 is kept, a is the 
universal parameter A = am and M R is the scale where the RH neutrinos are decoupled. 




1 



(Vy 2 3(^)33-/, 2 3 (3 + a 2 Klog 



Mgut_ 

Mr 



(2.9) 



8tt 2 



B. RGEs running 



The RGEs used by us are given in the Appendix A. We have paid much effort to keep the 
RGEs for the soft terms of lepton sector, which is relevant to our calculations, as complete 
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as possible. Specially speaking, both the diagonal and non-diagonal terms are kept in the 
RGEs. However, only diagonal terms of the scalar quark soft terms are kept because these 
non-diagonal terms are small and not relevant to our calculations. For the Yukawa sector, 
only the (3,3) elements of the diagonalized Yukawa matrices f t , ft,, f T and f U3 are kept. The 
evolution of the mixing matrix Vd from Mqut to Mr is ignored. 

The integration procedure consists of iterative runnings of the RGEs from the GUT scale 
to the low energy scale Mz and back for every set of inputs of m , mi , A Q , tan (5 and the sign 
of //, which are the universal scalar masses, gaugino masses, scalar trilinear parameter, and 
the standard vacuum expectation values ratio of the two Higgs fields and the Higgsino mass 
parameter, until the low energy gauge couplings and the Yukawa couplings are all correct 
within a given range. The parameters [i and B are given by low energy spontaneous breaking 
conditions Hl9lj20 |. 



The RGEs at Appendix A is given in the basis where f v is diagonal. The RH neutrinos 
are decoupled at Mr. Below Mr the basis is rotated to the basis where fi is diagonal. The 
MSSM RGEs are obtained by simply dropping the terms including f v and setting Vd equal 
to 1. The basis rotation leads to 

{m\) he iow = V D {m\) up Vb 

(A E ) below = V D (A E )V D . (2.10) 



Below Msusy — 250GeV r the SM beta functions are used 21 1. Threshold effects are taking 



into account by decoupling the corresponding particles at its running mass Q = m(Q). 



5m L 



Taking tan/3 = 2 ~ 10, we show the numerical results as follows 
/ 

(0.92 ~ 2.87) x 10~ 3 (-0.77 2.06) x 10 

(0.92 ~ 2.87) x 10~ 3 (0.97 ~ 3.0) x 10" 2 

^(-0.77 2.06) x 10^ 3 (0.97- 3.0) x 10- 2 



(3 + a 2 )m 2 Q (2.11) 



1 0.7 (2.42 ~ 8.12) x 1(T 4 

(2.42 ~ 8.12) x 10" 4 0.7 
I (-2.41 5.88) x 10~ 4 (2.5 ~ 8.67) 



x 10' 



where smaller tan /3 may give larger contribution. 



-2.41 ~ -5.88) x 10~ 4 
(2.5 ~ 8.67) x 10- 3 
0.7 



C. Low energy supersymmtric spectrum 



(2.12) 



At low energy the supersymmetric particle masses and mixing angles are obtained by 
diagonalizing the corresponding chargino, neutralino, scalar neutrino and scalar lepton mass 
matrices numerically. The slepton mass matrix is given by a 6 x 6 matrix as 



m LL m LR 1 
\ rrf LR m RR J 



where ttill, ^lr, wrr are all 3 x 3 matrices given by 



[m 



Lllj + (mj )6ij + M z(~2 + sin cos2(35ij , 



(m LR )ij = -((A E ) i:i + fitan/35 ij )mj 
{m RR )ij 



m R + mf — M| sin 2 W cos 2(3)6^ . 



(2.13) 



(2.14) 



The m 2 ^ is diagonal since only /; enters its RGE. The full sneutrino mass matrix has a 



12 x 12 structure. However, according to |T0| , p4]] if we only keep the first order of these terms 
perturbed by RH neutrino masses we can have a very simple structure, which is relevant to 
generation mixing, 



ml = m\ + \m\ cos 2(3 . 



(2.15) 



The mass matrices of chargino and neutrino are standard and given at Appendix B. 



III. ANALYTIC FORMULAS 



Fig. 2 gives the one loop diagrams relevant to the process r — ► /ry. The amplitude of 
this process can be written in the general form 



M = em i u j (p 2 )ia fll/ p'^(A L P L + A R P R )u i (p 1 )e fl (p 3 ) , 



(3.1) 



where Pl,r = |(1 =F 7s) are the chirality projection operators, i, j represent initial and 
final lepton flavor. The most convenient way to calculate A L and A R is to pick up the 
one loop momentum integral contributions which are proportional to Uj(p 2 )PL,RUi(pi)2pi • e 
respectively. The neutralino exchanging contribution is 



A™ = 



1 



1 



327r 2V v/2cosr m 



1 - 6k + 3k 2 + 2k 3 - 6k 2 log A; 



6(1 -ky 



A 



(n) 
R 



+ m X° a B jaa* A iaa I ~ k 2 + 2k log k 

mi (1 — k) 3 

l(n) 



AT (B^A), 



(3.2) 
(3.3) 



where k = m 2 /m 2 . A and B are the lepton-slepton-neutralino coupling vertices given in 

Xa 

Appendix B. The corresponding contribution coming from exchanging charginos are 

2 + 3k - 6k 2 + k 3 + 6k log k 



4 C) 



a 2 1 

j2 ryia* rvja 

32^ Z » Z » 



ml 



2a 2a 2M" 2 V cos 2 (3 " 



6(1 - k) 4 



m 



Xa 



z + z~ ' J 

V2M w cos(3 la 2a 



rrii 3 - 4k + k 2 + 2 log k 



mi 



'A-kf 



4 C) = 



92 gia* gja 



32tt 2 " 
m - 



V 1 

m u 

V<x L 



+ + *2 + 3k - 6k 2 + k 3 + 6k log k 



6(1 - k) 4 



V2M W cos (3 



3 - 4k + k 2 + 2 log k 
(1-A;) 3 



+ 



(3.4) 



(3.5) 



where k = m 2 -/m 2 . Mixing matrices Z„, Z + and Z are given in Appendix B. 

Xa a 



The branching ratio is given by 
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a 



BR(r - ni) = ^mtW + \A R \')/T T 



(3.6) 



where T T = 2.265 x W 12 GeV. 

Fig. 3 gives the diagrams contributing to the Z — > r/x process. Neglecting masses of the 
final fermions the amplitude is given as M = Eja^Mj, where Mj are 



M 1>2 = u(p 2 )YP L:R v( Pl )^ , 
M 3j4 = u(p 2 )Pl,rv(pi)pi ■ e 



(3.7) 
(3.8) 



respectively. The corresponding analytic expressions of are given in appendix C. 

The branching ratio for Z — > r/x is given by Br(Z — > /xr) = r(Z — > r ± /x T )/r^, where 
T z = 2.49(^6^ and 



r(Z - /xr) 



48ttM z 
1 

48ttM z 



(|a 1 | 2 + |a 2 | 2 )-2M| + (|a 3 | 2 + |a 4 | 



Ml 



(3.9) 



where £ represents the sum of /x, r spins and Z polarizations. 



IV. A SO(10) GRAND UNIFIED MODEL 

To give definite predictions we will work within a specific model. Before we turn to 
introduce this model it is worth noting that our calculations are not very model sensitive. 
That is the reason why we concentrate on the 2-3 generation mixing. The mixing between 
the first two generations may be very sensitive to different models and we will discuss the 
processes separately. From Fig.l, we know that r — /x mixing is mainly determined by the 
factor 5rhl 3 /ml, which, from Eq. ( pUD , depends only on the parameters (Vd)23 • (^0)33, 
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f U3 , M GUT and M R . M GUT k2x 10 b GeV is required by any supersymmetric grand unified 
model. In the context of SO(10) grand unified models /„ 3 is related to the top Yukawa 
coupling f t at the GUT scale. In addition, with the SuperK data of m Vr pa ^eV M R pa 
2 x 10 14 GeV A is roughly determined by the see-saw mechanism. Thus the only model sensitive 
parameter in fermion textures is (^0)23 ( (^0)33 is determined by unitary condition of the 
Vd matrix ). In most published SO(10) models the large — v T mixing is not produced 
dominantly by Mr mass matrix and thus produce large \i — r mixing, for example, we have 
(^0)23 • (^0)33 ~ (0.5 — 0.3) corresponding to #23 ~ (20° — 70°). Thus our discussions on 
these branching ratios are useful for estimating the predictions for such a class of unified 
models, although they may be different in details. 



We did our calculations within the model given by Albright et al ||15| . The model gives 
excellent predictions of quark and lepton masses and naturally explains the largeness of 

— v T mixing and the smallness of V c b- At Mgut, after the SO(10) breaking to the MSSM, 
the fermion mass textures are given by 



U 



7] 







-6/3 



\ 



N 



e/3 1 

r] 0^ 
e 
-e 1 j 



Mr, 



Mr. 



( 



D 





8 



8 




-e/3 



S'e^ a + e/3 1 

5 ^ 
5 a + e 
tie** -e 1 j 



M, 



D 



M. 



D 



(4.1) 



(4.2) 



where U, D, N and L are the up quark, down quark, neutrino and lepton mass matrices 
respectively. The most remarkable feature of the textures is the lopsided parameter a ~ 1 
at L and D. According to the SU(5) relation D = L T , the large — v T mixing due to 
the a term in L is related to the right-handed down quark mixing, which has no observable 
physical effects. The smallness of quark mixing V c b is determined by the parameter e ~ 0.1, 



11 



which translates to the right-handed lepton mixing. 

The model predicts the lepton sector 2-3 mixing with #23 ~ 63°, which leads to the only 
model structure sensitive quantity for our processes (Vb)23 • (^0)33 ~ 0.4. 

Taking all the fermion masses and Vkm elements at Mz given in Ref. as inputs, 
we calculate the corresponding values at Mqut with several values of tan (3 by integrating 
the RGEs and fit the parameters in Eq. (4.1) and (4.2). We find these parameters are not 
sensitive to the supersymmetric parameters, except that Mu becomes larger as taking small 
tan (3. To keep the predicted neutrino masses in the correct range we take Mr = 5 x 10 14 GeV 
when tan (3 = 2 and Mr = 2 x 10 1A GeV when tan (3 = 5, 10 in our later calculations. 



V. NUMERICAL RESULT AND DISCUSSION 

The relevant parameters on predicting the branching ratios are the universal supersym- 
metry soft breaking parameters, m , mi, A , tan (3 and the lepton sector mixing matrix Vd 
defined in Eq. (2.3). We will present dependence of the branching ratio of the processes 
r — > /i7 and Z — > t/i on the soft parameters in the model of Jl5|]. In our calculations the soft 
parameters are constrained by various theoretical and experimental considerations |23[ , such 
as the LSP of the model should be a neutralino and masses of all supersymmetric particles 
be beyond present mass limits and so on. 

In Figs. 4~7, we plot the branching ratio of r — > /ry as functions of tuq for several 
sets of other soft parameters. Two experimental bounds are plotted in every figure, which 
correspond to the present experimental limit BR(t — > fi^) < 1.1 x 10~ 6 [24] and the 



expected sensitivity of future experiments BR{t — > /17) < 1.0 x 10~ 9 |T(|. The general 
trend of BR(t — > ^7) is decreasing dramatically as m increases. In Fig. 4, the branching 
ratio is plotted for mi = lOOGe^, tan/5 = 5 and A = 2m , m , 0, taking /i either positive 



12 



or negative. That branching ratio decreases with a = Aq/ttiq is easily understood due to the 
Eq. (1). It is interesting to note that although this process can not be observed at present, 
it will be detectable in the future experiments in a large part of the soft parameter space. 

Fig. 5 presents the branching ratio for mi = 100, 200, AOOGeV and A Q = m , tan (3 = 5. 
We can find that the branching ratio is also very sensitive to the parameter mi. It will 
decrease quickly with increasing mi . If mi > 400GeV, this process will be non-observable 

2 2 

as displayed in the figure. In Fig. 6, we plot the branching ratio for mi = lOOGeV^, 
tan/3 = 10 and A = ±2m , ±m , 0. We find that the sign of A Q is not very important 
in the order of magnitude estimate. However, the sign of Aq is still relevant to the precise 
predictions. Fig. 7 plots the branching ratio at different tan (3. The branching ratio increases 



as tan (3 becomes large. It is explained in Ref. |1^] that the dominant term of the amplitude 
for r — ► /i7 is proportional to tan f3. We note that at tan (3 = 2 the sign of \x can be 
significant. The branching ratio when fi is negative is much smaller than that when fi is 
positive. 

In Figs. 8~10, we present the quantitative results of the branching ratio of the decay 
Z — > TjJ,. We find in most parameter space this process can not be detected by the Giga-Z 
option [17]. However, we can still find several interesting features in this process which are 



different from the r radiative decay process. The most remarkable feature about this LFV 
process is that its branching ratio becomes large at first and then approaches to a constant 
when mo is increasing, which gives a sharp contrast with the process r — > /ry. Thus this 
process may have advantage over the r — > /ry process in some region of parameter space. 
The reason for this different behavior between the two processes can be traced back to the 
different coupling structures between the processes as shown in Eqs. (3.1) and (3.7), (3.8). 
The magnetic structure of r decay determines that its amplitude is inversely proportional 
to the sfermion masses square, whereas there is a vector current coupling in Eq. (3.7) which 
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determines the Z — > /it amplitude's trend as increasing tjiq. 

Fig. 8 shows the branching ratio as a function of m for mi = lOOGeV, tan/5 = 5 
and A = 2m , m , with both positive and negative sign of /i. The sign of \x is quite 
irrelevant when m > 500GeV ". Fig. 9 displays the same function for different mi values, 
100, 200, 400GeV, with = 2mo and tan/3 = 5. The dependence of BRiZ — > r/i) on a 
and mi is similar to that of BRir — ► ^7). We note again that the sign of /i is irrelevant when 

mi > 200GeV. Fig. 10 gives the branching ratio for tan/? = 2, 5, 10 with A = 2.5m , 

2 

mi = 150GeV. We can see another important feature of this LFV process that it is more 

2 

favorable for the small tan (3 value in contrast with the r — > /J7 process. At the extreme case, 
the branching ratio may access 1 x 10~ 8 , which may be detectable. The relationship between 
the branching ratio and tan (3 is easily understood. For small tan (3 large non-universal soft 
mass term <5m| 3 will be produced due to a large f V3 in Eq. ( |2.9| ). We note that when tan/3 
becomes smaller, the BRiZ — > r/i) increases quickly. 

In summary, according to our calculations we find that the r — ► is more feasible than 
Z — >• r/i to reveal charged lepton flavor violation in the context of SO(10) SUSY-GUTs. 
In most parameter space r — ► /ry has a branching ratio that can be detected in the future 
experiments JTIJ . The Z — > r/i is hopeful to be detected in a small parameter space. We note 
the remarkable feature of the Z decay process that its dependence on the supersymmetry 
soft parameters tan (3 and mo is opposite to that of r decay. Thus it can supplement the 
LFV search besides the r — > /ry process. 



VI. SUMMARY AND CONCLUSION 

In this work, we present the dependence of the branching ratios of two charged lepton fla- 
vor violation processes r — > 7/7 and Z — ► r/i on the supersymmetry soft breaking parameters 
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in the context of SO (10) grand unified models with the "lopsided" texture of mass matrices 
for charged leptons and down quarks. We expect these processes may be detected in the 
future experiments. The first process is more hopeful to be observed. The second process 
may offer useful information about the soft parameters if it is also observed. The different 
behaviors of the two processes depending on tan/? and mo implies that the simultaneous 
study of the two processes will be interesting. 

We expect the study of charged lepton flavor process may provide another window of 
high energy physics besides the neutrino oscillation study. The combined study of the LFV 
processes and neutrino oscillation may feed light on the sector of right-handed neutrinos, 
which may be necessary in a model to naturally explain the small neutrino masses. 
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Appendix A 



In this appendix we give the renormalization group equations of the MSSM suplemented 



with RH neutrinos |p~2| , [T3|] . The two-loop RGEs of the gauge couplings can be found in 
many literatures which will not be affected by the presence of the RH neutrinos. We give 
one-loop RGEs of the Yakawa coupling matrices and the soft terms which are affected by 
the presence of RH neutinos. In this appendix we denote the Yukawa couplings of up quark, 
down quark, lepton and neutrino as U, D, E and iV respectively. Denote the soft terms as 
Ajj ■ U = U a and so on. The RGEs of the Yukawa couplings are 

1 



dU 

~dt ~ 16tt 2 



-E Ci c/ 2 + WU ] + DD^ + Tr[3UU^ + NN^} 



U 



(A.l) 
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dD 

lit 
dE 

~dt 
dN 

~dt 



16tt 2 
1 

16tt 2 
1 

16/T 2 



-£c^ 2 + 3DD ] + UU f + Tr[3DD ] + EE^] 
-Zc'-g 2 + 3E-E t + AW* + Tr[3DD^ + EE^] 
-Ec^'ft? + 3iViV f + EE 1 " + Tr[3f/f/ f + iViV f ] 



E 
N 



(A.2) 
(A.3) 
(A.4) 



where t = logQ, c, = (if, 3, f ), c, = 3, f ), cj = (|, 3, 0), = (|, 3, 0). The RGEs of pt 
and soft parameters of Higgs sector are 
djj, 1 



dt 16tt 2 
dB 2 



-Eq + Tr[3UU ] + 3/^ + + NN^\ 



dm 2 rr 



dt 16tt 2 
2 



dm 2 HD 



16tt 2 



-Ec'-'g^Mi + Tr[3UU A + 3DD A + ££ A + NN A ] 



^g 2 M 2 + 3Tr[U(M* L + M* R )rf + m^C/C/t + ^[/t] 



dt 



+Tr [Mq l NN^ + NM~ R N^ + m 2 Hu NN^ + A^A^] 



-^g 2 M 2 + 3rr[£>(JM3 L + M^J^t + m 2 HD Drf + D^] 



16tt 2 



+Tr[E(M 2 L + Ml)E^ + m 2 Ho EE^ + E A E A ]\ . 



(A.5) 
(A.6) 

(A.7) 
(A.8) 



Mi in the above expressions are the gaugino masses whose RGEs are same as those in the 
MSSM. 



Then we give RGEs of soft mass terms of lepton sector 

< ^§ ±= \L 2 + l^ NNtM l + M 2 L NN^] + \\EE^M\ + M\EE^\ 

+EM\E ] + m 2 HD EE^ + E A E A + NM 2 R N^ + m 2 Hu NN^ + N A N A ] , 



d_M\ 
dt 

dt 



167T 2 

2 

167T 2 



-^g 2 M 2 + E^EM\ + M\E^E + 2(£ t Mf £ + m 2 Ho E^E + E A E A ) 
o 

N^N + N A N A ) 



(A.9) 
(A.10) 
(A.ll) 



The RGEs of trilinear terms of lepton sector are 
dA E 1 



dt 



16tt 2 



2Zc'-g 2 M i + 2Tr{3A D DD ] + A E EE ] ) + 2A N NN* 
16 



dA 



N 



dt 



+(5££ f + NN r )A E + A E (EE* - NN r )] , 

^ \-TLclg\Mi + 2Tr(3AuUtf + A N NN*) + 2A E EE* 

+ (5A^V f + EE r )A N + A N {N^ - EE r )] . 



(A.12) 



(A.13) 



In the basis where N is diagonal and only keep the third family Yukawa coupling constants 
ft, fb, fv 3 and f T , then the RGEs are simplified as 
1 



djt 
dt 
dfo 
dt 

dt 

#r 

dt 
djj, 

~dt 
dB 

~dt 

dm 2 Hu 

dt 



dm 2 Ho 
dt 



'Ml 
dt 



( dm R 
\ dt , 

' drriv^ 



dt 



dA 

dt 
dA b 

~dT 



+ QfZ + ft + fr\ fb , 

-xc^+ifi+vt+m 



-Zc:g 2 + 4f T +3fZ + f\V 33 \ 2 



16tt 2 
1 

16^2 
1 

16V 2 
1 

167T 2 
1 

16tt 2 
2 

167T 2 " 

2 

167T 2 



-Scr^ + 3/ 2 + 3/ 2 + / 2 + 



fl>3 1 
fr , 
fi 



EcJ^M, + 3/ 2 A + 3ftA b + £4, + / 2 A T 



-Xc: 9 *M? + 3f?(Mi+M- 



m 



+ A 2 ) 



i/ 3 ; 



-E C ;"^ 2 M 2 + 3/ 2 (M? + Mr + m 2 



+ f 2 T (M 2 fL + Ml+mj lD +Al) 
2 



A 2 ) 



16tt 2 



Sc, <? 2 M 2 )^ + VZVvft ( -((M|)jj + (MQv) 



^3 H^f),^/) K/,,^),,/);, K .Vj>, rW>, ; + (.4 (/ ),: { (^). /:5 



167T 2 

2 

167T 2 

2 

ilv 2 

2 

16V 2 
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5 ^ 2 M 2 + 2/ 2 ((M|) ll + (M|)jj + m 2 ^ + A 2 )fe 



2/ 2 3 ((Ml) ii + (M?> + m 2 ¥t/ +A 2 ) 



-Ec^M, + 6/ 2 A + f h A h + / 2 A, 3 



-Ec^M, + 6fUb + /X + / t % 
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Ji3 



(A.14) 
(A.15) 
(A.16) 
(A.17) 
(A.18) 
(A.19) 

(A.20) 

(A.21) 

(A E ) u (A E )jj 

, (A.22) 
(A.23) 

(A.24) 

(A.25) 
(A.26) 



dA E 



' dA v " 



16vr 2 



;-2S c ;'^ 2 M l + Qf 2 b A b + 2f T A T )5 l3 + \];y :ij f 2 (MA ,:),,, + {A E ) U )+ 



f* 3 (2(A u ) i3 S j3 +(A E ) 3j 6 i3 -(A E ) i3 d j3 )\ , (A.27) 
j -2E C ;"^ 2 M i + 6f?A t + 2/ 2 3 A 3 )% + V^y 3 i/r(2(A B )« + - (A„)«)+ 



167T 2 



flMA v ) 3j 8 l3 + {A v ) l3 5 ]3 ) 



(A.28) 



The matrix in the above equations refers to Vp defined in Eq. (2.3). Below Mr, RH 
neutrinos are decoupled and the RGEs of the MSSM are used. This can be achieved by 
setting f U3 = and V!y = Sij in the above expressions. 

Appendix B 

In this appendix we list the relevant pieces of the Lagrangian and conventions which we take 
in our calculations fll9|. The Lagrangian pieces are 



2 cos &w 

C Tlz = (\ZHZT - sin 2 My) (5 K l 3 )Z 

cos a y/ \l / 



r 9% n u 

92 Xh» 



[Z+*Z+P R + Z^Z U *P L + Sa cos26 w ] X - Z„ , 



2 cos 6 W ' 



K^N ~ ^Nl^L — N ~ Z 'N Z 'N ) r R 

~C ij P L + D ij P R 



Y° Z 
Xj M 7 



zf*P R - 



mii 



V2M W cos /3 

zL\zlj + 



ZI J X-KJ 



+ h.c. 



(B.l) 
(B.2) 
(B.3) 

(B.4) 
(B.5) 



Z 2 J cot 9 W ) - cot 9 W ™ 11 - Zf y Z z A P L 
M w cos p J 



- Uz^zf + cote w -^z»zf) p R ] i 1 t* + h.c. 

\ M w cos p ) 



y2 cos 8w 



'A^Pl + S I<J 'P fl ] I 1 !* + h.c 



(B.6) 



The abbreviations defined in the above expressions will be used in the next Appendix. 
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The mixing matrices Z in the above expressions are given below. The scalar lepton and 
scalar neutrino mass matrices are given in Sec II. The corresponding mixing matrices are 
defined as 



Z\vrc~Z L = diag (m~J , 
Z\m\Z v = diag (m|.) , 



i = 1...6 

i = 1, 2, 3 



(B.7) 
(B.8) 



The mass matrix of charginos is 



The mixing matrices Z^ is defined as 



m 2 \/2Mw sin (3 

\plM w cos (3 /I 



[Z ) T M X Z + = diag (m Xl ,m x 2 ) 



The mass matrix for neutralinos is 

mi 

m 2 
—Mz cos (3 sin 9w M z cos ft cos 9w 
M z sin (3 sin 0^ — M z sin /3 cos 6^ 



M x o = 



—M z cos /3 sin 6 ] w 
Mz cos /3 cos 9w 




(B.9) 



(B.10) 



M z sin /5 sin 9w 
—Mz sin /3 cos 6w 
-fi 




(B.ll) 



which is diagonalized by 



ZJ^M x oZn = diag (jn x o, m x o, m x o, m x o^j 



;b.i2) 



Appendix C 



In this appendix we give the analytic expressions of defined in Sec III. 
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a\ is given by 



a-y 



92 92 

2 cos 9w 16-7T 2 



(ai(a) + ai(fe) + ai(c) + ai(d) + 2ai(e) + 2ai(/)) , 



(C.l) 



where ai(i) is coming from the coresponding Feynman diagram. They are given by 



ai(a) 
a>i(b) 
ai(c) 

ai(d) 
ai(e) 



= |Z 1 + J 2 ^Zr*(-2)C o 



Z+ * 7+ via via 11 



+ 5 ab cos29 w )(B - 2C 00 



(C.2) 
(C.3) 

{{Z- h Z^ + 5 a6 cos2^)m x -m x - - (Z+*Z+ + S ab cos29 w )ml a ) C 

(C.4) 



- tan 2 9 W (^Z^Z 1 / - sin 2 9 w 5 a ^ (-2)A iaa * A^ a C m , 



-- tan 2 e w A iaa *A jab [2D ab C 00 + (C ab m x om x » - D ab m 2 L ) C - D ab B 
-2(0.5 - sin 2 e w )\Z+\ 2 Z?Zi a *(B + B,) 
— tan 



! 9 W (0.5 - sin 2 9 w )A iaa *A jaa {B + B x ) , 



(C.5) 
(C.6) 
(C.7) 



where i, a and a represent the flavors of lepton, slepton or sneutrino and chargino or neu- 
tralino respectively. a 2 is given by 



a 2 



92 92 

2 cos 9 W 167T 2 



(02(b) + a 2 (d) + 2a 2 (/)) , 



(C.8) 



where 



a 2 (6) = -t^e w qzl a Zi P - sin 2 9 w 5^)(-2)B^B^ a C 00 , 
a 2 (d) = -itaQ0 w S iao *B J '°* [2C a6 C o + (/) ab m x om x o - C afe m?J C - C ab B 
02(f) = tan 2 9 W sin 2 9 w B iaa * B^ aa (B + fix) , 



(C.9) 
(CIO) 
(C.11) 



Then we have 



92 



9l 



2 COS#jy 167T 2 



(oi(6) +Oi(d)), i = 3, 4 



(C.12) 
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and 



03(b) = 2 tan 2 Owi-Z^Z 1 / - sin 2 6 w 5 a P)B iaa * A j/3a m x «(C + C x + C 2 ) , (C.13) 

a 3 (d) = tan 2 B w B iaa * A ja \m x »C ah C x + m x oL> a6 C 2 ) , (C.14) 

a A (b) = 2 tan 2 d w (-Z^Z^ - sin 2 ^<?^)A^*B^ a m xg (C r + C x + C 2 ) , (C.15) 

a 4 (d) = tan 2 w ^ iaa *5 ia6 (m x o J D a6 C 1 + m x oC a6 C 2 ) . (C.16) 



All the coupling vertices A, 5, C, D are defined in the Appendix B. The I?o,i and 
Co,i,2,oo are the standard two-point and three-point functions with its defination given in 
the program LoopTools ||25|| . The arguments of function B from Fig. (c), (d), (e), (f) are 
(Mf, m 2 _, m 2 _), (Mf, m 2 , m 2 ), (0, mj , m 2 -) and (0, m~ , m 2 ) respectively. The 

Xa Xf, Xa Xjj a Xa '■a Xa 

arguments of functions C from Figs, (a), (b), (c), (d) are (0, Mf, 0, m 2 _, m l a , m 1 a )i 
(0, Mf, 0, m 2 , mf , mf ), (0, M|, 0, m? Q , m 2 , m 2 ), (0, Mf, 0, mf , m 2 , m 2 ) 
respectively, where the external fermion masses have been set to zero. 
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FIGURE CAPTIONS 



FIG. 1 Feynman diagrams contributing to r — > pq. The universal scalar neutrino masses 
become non-degenerate and mixed by the non-diagonal Yukawa couplings. This effect 
gives low energy scalar neutrino mixing by RGEs running. 

FIG. 2 The Feynman diagrams contributing to the process r — > /ry. 

FIG. 3 The Feynman diagrams contributing to the process Z — > /xt. The other two self- 
energy diagrams coming from r legs are omitted. 

FIG. 4 Branching ratio of r — > /17 as a function of itlq for mi = lOOGeV, tan/3 = 5 and 
A = 2m , mo, 0. The solid line is for /x > and the dashed line is for /x < 0. 

FIG. 5 Branching ratio of r — > /17 as a function of m for A = m , tan/3 = 5 and 

mi = lOOGeV, 200GeV, 400GeV. The solid line is for a > and the dashed line is 

2 

for /x < 0. 

FIG. 6 Branching ratio of r — > ^7 as a function of m for mi = lOOGeV, tan/3 = 10, 
/i > and A = ±2m , ±m , 0. The solid line is for A > and the dashed line is 
for A < 0. 

FIG. 7 Branching ratio of r — ■> 7x7 as a function of m for A = m , mi = lOOGeV, and 
tan/3 = 2, 5, 10. The solid line is for /i > and the dashed line is for /x < 0. 

FIG. 8 Branching ratio of Z — > /xr as a function of m for mi = lOOGeV, tan/3 = 5 and 
A = 2m , m , 0. The solid line is for /x > and the dashed line is for /x < 0. 

FIG. 9 Branching ratio of Z — > /xr as a function of m for A = 2m , tan/3 = 5 and 

mi = lOOGeV, 200GeV, 400GeV. The solid line is for a > and the dashed line is 

2 

for /x < 0. 
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FIG. 10 Branching ratio of Z — > fir as a function of m for A = 2.5m , mi = lSOGeV, 
and tan (3 = 2, 5, 10. 
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